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We present measurements of branching fractions and charge asymmetries for six B-meson decay 
modes with an r/ or rf meson in the final state. The data sample corresponds to 232 million BB pairs 
collected with the BABAR detector at the PEP-II asymmetric-energy e + e _ B Factory at SLAC. We 
measure the branching fractions (in units of 10 -6 ): B(B + —> r/n + ) = 5.1±0.6±0.3, B(B + —pK + ) = 
3.3±0.6±0.3, B(B° pK°) = 1.5±0.7±0.1 (< 2.5 at 90% C.L.), B(B+ -> pp+) = 8.4±1.9±1.1, 
B(B° — > r/ui) = 1.0 ± 0.5 ± 0.2 (< 1.9 at 90% C.L.), and B(B + —^ ?/7r + ) = 4.0 ± 0.8 ± 0.4, where the 
first uncertainty is statistical and second systematic. We also determine the charge asymmetries for 
the charged modes: A c h{B + —> r/n + ) = —0.13±0.12±0.01, A c h{B + —> r/K + ) = —0.20±0.15±0.01, 
Ach(B+ -> pp+) = 0.02 ± 0.18 ± 0.02, and A c h(B+ -> p'n + ) = 0.14 ± 0.16 ± 0.01. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 


Charmless B decays are becoming increasingly useful 
to test the accuracy of theoretical predictions, for exam¬ 
ple based on QCD factorization 1, [2j or flavor SU(3) 
symmetry 0, [ij. In this Letter we present measurements 
of branching fractions and, when applicable, charge 
asymmetries for six charmless B decays: B + —> pn + , 
B + —> rjK + , B° —> pK°, B° —> put, B + —> r/p + and 
B + —> 7 / 7 r + , of which the last four were not observed 
before |5, [fj, 7j, |8]. Some of these decays may proceed 
through CKM-suppressed b —> u and loop (“penguin”) 
b —> s transitions with amplitudes of comparable size. 
Interference between these amplitudes can lead to direct 
CP violation measurable in charge asymmetries A c h Q- 
The measured branching fractions and charge asymme¬ 
tries may also be sensitive to the effect of non-Standard- 
model heavy particles entering the loop 0 . 

Charmless B decays with kaons are usually expected 
to be dominated by b —> s loop amplitudes, while b —> u 
tree amplitudes are typically larger for the decays with 
pions and p mesons. However, the B —► pK decays are 
especially interesting since they are suppressed relative 
to the abundant B —> p'K decays due to destructive 
interference between two penguin amplitudes E3- The 
CKM-suppressed b —> u tree amplitudes may interfere 
significantly with b —> s penguin amplitudes of simi¬ 
lar sizes, possibly leading to large direct CP violation 
in B + —> ?y p + , B + —> ? 77 r + and B + —> p l n + 11]; nu¬ 
merical estimates are available in a few cases Ulna. 
We search for such direct CP violation by measuring the 
charge asymmetry A c h = (r _ — T + )/(r _ + T + ) in the 
rates T 1 * 1 = T(I3 ± —> f ± ) for each charged final state f ± . 

Finally, phenomenological fits to the branching frac¬ 
tions and charge asymmetries of charmless B decays can 
be used to understand the relative importance of tree and 
penguin contributions and may provide sensitivity to the 
CKM angle 7 00 

The results presented here are obtained from extended 
unbinned maximum likelihood (ML) fits to data collected 
with the BABAR detector [H at the PEP-II asymmetric 
e + e _ collider 0 located at the Stanford Linear Accel¬ 


erator Center. The analysis uses an integrated luminos¬ 
ity of 211 fb -1 , corresponding to 232 million BB pairs, 
recorded at the T(4S') resonance (center-of-mass energy 
7 ^= 10.58 GeV). 

Charged particles are detected and their momenta 
measured by a combination of a vertex tracker consisting 
of five layers of double-sided silicon microstrip detectors 
and a 40-layer central drift chamber, both operating in 
the 1.5-T magnetic field of a superconducting solenoid. 
We identify photons and electrons using a CsI(Tl) elec¬ 
tromagnetic calorimeter (EMC). Charged particle identi¬ 
fication (PID) is provided by an internally reflecting ring 
imaging Cherenkov detector (DIRC) covering the central 
region of the detector, the average energy loss ( dE/dx) in 
the tracking devices and by the EMC. A K/tt separation 
better than two standard deviations (<r) is achieved for 
all momenta. 

We select 77 , 77 ', u>, K ° and 7 r° candidates through the 
decays 77 -> 77 ( 77 -^), 77 -> 7T + 7r _ 7r° (? 7 37r ), 77 ' -> 77 77 7t + 7t _ 
K™). o' -> P °7 Wp-y), U -> TT+TT-TT 0 , K° -> TT+TT" and 
7 r° —> 77 . We impose the following requirements on the 
invariant mass in MeV of the particle candidates’ final 
states: 490 < m 77 < 600 for 77 77 , 520 < < 570 for 

773 ^, 910 < (nirfKir, rripj ) < 1000 for 77 ', 735 < m „ v < 825 
for oj, 510 < m n7r < 1070 for p°, 470 < m n7r < 1070 
for p + , 486 < < 510 for K° and 120 < TTl-yry < 

150 for 7 r°. These cuts are loose for the invariant mass 
variables used in the ML fit, and tight for those that 
are not. For K® candidates we require at least 3er three- 
dimensional separation between the decay vertex and the 
e + e _ collision point. For the vector resonances to and 
p + we also use the hclicity-frame decay angle 0h- The 
helicity frame is defined as the vector-meson rest frame 
with polar axis along the direction of the boost from the 
B rest frame. For u>, Oh is the polar angle of the normal 
to the decay plane, and for p it is the polar angle of the 
charged daughter momentum. We define H = cos Oh and 
require —0.75 < H < 0.95 for p + . 

All tracks from resonance candidates are required to 
have PID consistent with pions. For the B + decays to 
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? 77 r + , r]K + and ? 7 7 7 r + , the primary charged track must 
have an associated DIRC Cherenkov angle within 3.5cr of 
the expected value for either a it or K hypothesis. The 
discrimination between primary 7 r and K is performed in 
the ML fits. 

A B-meson candidate is characterized kinematically 
by the energy-substituted mass toes = (i s ~ Ps ) 5 and 
energy difference A E = Eb — \yfs, where ( Eb,Pb ) is 
the B-meson 4-momentum vector, and all values are ex¬ 
pressed in the T(4S) frame. Signal events peak at zero 
for A E, and at the B nominal mass for toes- The res¬ 
olution on A E (toes) is about 30 MeV (3.0 MeV). We 
require \ AE\ < 0.2 GeV and 5.25 < toes < 5.29 GeV. 

Backgrounds arise primarily from random combina¬ 
tions in continuum e + e _ —> qq (q = u,d,s,c) events. 
To reject these events we make use of the angle 9t be¬ 
tween the thrust axis of the B candidate in the Y(AS) 
frame and that of the rest of the charged tracks and neu¬ 
tral clusters in the event. The distribution of | cos0t| is 
sharply peaked near 1 for combinations drawn from jet¬ 
like qq pairs, and nearly uniform for the almost isotropic 
B-meson decays; we require |cos$t| < 0.9 (< 0.65 for 
if pl 'K + ). Further discrimination from continuum in the 
ML fit is obtained from a Fisher discriminant T that is 
described in detail elsewhere @- 

Where necessary, we use additional event selection cri¬ 
teria to reduce BB backgrounds from several charmless 
final states. Specifically, we require that photons have 
energies in ranges uncharacteristic of these backgrounds 
and, in ? 7 77 cl> and 77 77 A'°, we eliminate ?y 77 candidates 
that share a photon with any 7r° candidate having mo¬ 
mentum between 1.9 and 3.1 GeV in the Y(4S) frame. 

Multiple candidates are found in less than 30% of the 
events, in which case we choose the candidate with the 
smallest value of a y 2 constructed from the deviations of 
the daughter resonance masses from their nominal values. 

We use Monte Carlo (MC) simulation 0 for an initial 
estimate of the residual BB background and to identify 
the few (mostly charmless) decays that may survive the 
candidate selection and have characteristics similar to the 
signal. We find these contributions to be negligible for 
several of our modes. Where they are not negligible, 
namely for 77 77 7t + , ? 7 77 A' + , r\p + and q' pry n + , we include a 
component in the ML fit to account for them. 

We obtain yields and A c h for each decay chain from a 
ML fit with the following input observables: A E, Toes, 
T , and m res (the mass of the 77 , 77 ', p + , or ui candidate). 
For uj and p + decays we also use H and, for charged 
modes with a primary charged track, the PID variables 
Sit and Sk , defined as the number of standard deviations 
between the measured DIRC Cherenkov angle and that 
expected for pions and kaons, respectively. 

For each event i, hypothesis j (signal, continuum back¬ 
ground, BB background), and flavor k (primary 7 r + or 


AT + ), we define the probability density function (PDF) 

V) k = Vj{rriES i )'Pj(AEl [, ^])B,(^)B J (m( es [,W]) (1) 

The bracketed variables S and H pertain to modes with 
a primary charged track or vector resonance daughters, 
respectively. Known correlations between A E k and Sk, 
and between m Tes and 7Y, are included in the PDF. The 
likelihood function is 


N 

C = exp (- Yjk) n 
J,fc i 


j,k 


( 2 ) 


where Yjk is the yield of events of hypothesis j and fla¬ 
vor k, to be found by maximizing C. N is the number of 
events in the sample. Free parameters of the fit are the 
signal and background yields, qq background PDF pa¬ 
rameters (see below), and for charged modes the signal 
and qq background charge asymmetries. 

For the signal and BB background components we de¬ 
termine the PDF parameters from simulation. For back¬ 
ground from continuum (and non-peaking combinations 
from B decays) we obtain the PDF from (toes , A E) side¬ 
band data for each decay chain, before applying the fit 
to data in the signal region; we refine this PDF by let¬ 
ting as many of its parameters as feasible free to vary 
in the final fit. We parameterize each of the functions 
^sig(wEs)) 'Psig(ABfc), Vj[T), V{S k ) and the peaking 
components of Vj(rn res ) with either a Gaussian, the sum 
of two Gaussians or an asymmetric Gaussian function 
as required to describe the distribution. Slowly varying 
distributions (mass, energy or helicity-angle for combi¬ 
natorial background) are represented by one or a com¬ 
bination of linear, quadratic and phase-space motivated 
functions [t|. The peaking and combinatorial compo¬ 
nents of the u> and p + mass spectra each have their own 
H shapes. Control samples with similar topologies as our 
signal modes (e.g. B —> D(Ktttt)tt) are used to verify or 
adjust the simulated resolutions evaluated from MC. 

Before applying the fitting procedure to the data we 
subject it to several tests. In particular we evaluate pos¬ 
sible biases in the yields from our neglect of small resid¬ 
ual correlations among discriminating variables in the 
PDFs. This is achieved by fitting ensembles of simu¬ 
lated qq experiments drawn from the PDF into which we 
have embedded the expected number of signal and BB 
background events, randomly extracted from the fully 
simulated MC samples. The measured biases are listed 
in Table |U 

The branching fraction for each decay is obtained from 
the measured yield, corrected for the fit bias and for the 
selection efficiency, and the number of BB pairs. We as¬ 
sume equal decay rates of the T(4S) to B + B~ and B°B°. 
In Table |U we show for each decay mode the measured 
branching fraction together with the event yield and effi¬ 
ciency, and A c h when applicable. The purity is the ratio 
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TABLE I: Fitted signal yield Ys in events (ev.), estimated purity P, measured bias (see text), detection efficiency e, daughter 
branching fraction product (]”[ Bi), significance S (with systematic uncertainties included), measured branching fraction B, and 
signal charge asymmetry Ach for each mode. The quantities in parentheses are 90% C.L. upper limits. 


Mode 

Y s (ev.) 

P (%) Bias (ev.) 

e (%) n & (%) 5 (a) 

B (10“ 6 ) 

*A c h 

r)- yl 7r + 

153^22 

30 

+7 

33 

39 

7.9 

4-8+g® 

-0.04 + 0.14 

V3^ + 

76±j s 

32 

+6 

24 

23 

5.6 

5.6tli 

-0.32 + 0.20 

J77T + 






9.7 

5.1 ± 0.6 ± 0.3 

-0.13 ± 0.12 ± 0.01 

7]~ n K + 

nets 

29 

+8 

32 

39 

6.1 

3.6 ±0.7 

-0.19 + 0.16 

V3 nK+ 

qy+13 

-12 

24 

+5 

23 

23 

2.8 

2 6 +11 

-0.22 + 0.33 

r)K+ 






6.7 

3.3 ± 0.6 ± 0.3 

-0.20 ± 0.15 ± 0.01 

Tj-y-yK 0 

17 + ? 

27 

+3 

28 

14 

2.3 

1 fi +10 
^-o.g 


V3 nK 0 


28 

+1 

21 

8 

1.4 

1 1 +1 - 3 
^-L-o.g 


r/K 0 






2.6 

1.5 ± 0.7 ± 0.1 (< 2.5) 


TJ'y'yLO 

13 ±1 

32 

+1 

14 

35 

2.5 

l-ltois 


r]3nUJ 

2 tl 

6 

-1 

11 

20 

0.6 

0-6t^ 


r]LO 

126±32 





2.5 

1.0 ± 0.5 ± 0.2 (< 1.9) 


7 hi P + 

12 

+18 

16 

39 

3.7 

7.3j£| 

0.10 + 0.23 

r)3nP + 

65t£ 

15 

+3 

11 

23 

3.4 

10-6 til 

-0.14 + 0.31 

VP + 






4.7 

8.4 ± 1.9 ± 1.1 

0.02 ± 0.18 ± 0.02 


69±H 

42 

+9 

27 

18 

5.6 

5.5±i? 

0.09 + 0.18 

»?P7 7r+ 

30±i 6 5 

13 

+9 

17 

30 

1.4 

1.8l(J 

0.58 ± 0.44 

ri'n + 






5.4 

4.0 ± 0.8 ± 0.4 

0.14 ± 0.16 ± 0.01 


of the signal yield (Ys) to the effective background plus 
signal (FJ ff +Ys), which we estimate as the square of the 
uncertainty in the signal yield (YJ ff + Ys = 4J- 

The statistical uncertainties in the signal yield and A c h 
are taken as the change in the central value when the 
quantity —2 In £ increases by one unit from its minimum 
value. The significance is taken as the square root of the 
difference between the value of — 21n£ (with systematic 
uncertainties included) for zero signal and the value at 
its minimum. 

For each mode the measurements for separate daughter 
decays are combined by adding the values of — 21n£ as 
functions of branching fraction, taking proper account of 
the correlated and uncorrelated systematic uncertainties 
described below 0. For r]LO and r]K° we quote 90% confi¬ 
dence level (C.L.) upper limits, taken to be the branching 
fraction below which lies 90% of the total of the likelihood 
integral in the positive branching fraction region. 

In Fig. ^ we show projections onto ju-es and A E of 
subsamples enriched with a mode-dependent threshold 
requirement on the signal likelihood (computed without 
of the variable plotted) that optimizes the sensitivity. 

Most of the uncertainties arising from lack of knowl¬ 
edge of the PDFs have been included in the statistical 
uncertainty since most background parameters are free 
in the fit. For the signal the uncertainties in PDF pa¬ 
rameters are estimated from the consistency of fits to 
MC and data in control modes with similar final states. 
Varying the signal PDF parameters within these errors, 
we estimate the mode-dependent uncertainties due to the 
signal PDFs to be 1-8 events. We verify the validity of 



FIG. 1: The B candidate tues (left) and A E (right) projec¬ 
tions obtained with a cut on the signal likelihood (see text) 
for B + —> pp + (a, b), B + — > r\ 7r + (c, d), and combined 
B + —> 7]-k + and B + —> r/K + (e, f). Points with uncertainties 
represent the data, solid curves the full fit functions, dashed 
curves the background functions and the dotted curves the 
background plus signal r/K + functions. 

the fit procedure and PDF shapes by demonstrating that 
the likelihood for each fit is consistent with the distribu¬ 
tion found in simulation. 

The uncertainty in the fit bias correction is taken to 
be half of the correction itself. Similarly we estimate 
the uncertainty from modeling the BB backgrounds by 
taking half of the difference between the signal yield fitted 
with and without the BB background component. 
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Uncertainties in our knowledge of the reconstruction 
efficiency, found from auxiliary studies on inclusive con¬ 
trol samples 0, include 0.6% per primary track, 0.8% 
per track from a resonance, 1.5% per photon, and 2.1% 
for a K®. Our estimate of the systematic uncertainty in 
the number of BB pairs is 1.1%. Published data [16] pro¬ 
vide the uncertainties in the 73-daughter product branch¬ 
ing fractions (1-3%). The uncertainties in the efficiency 
of the event selection are 1% (4% in B + —> p' pl 'n + ) for 
the requirement on cos(9t and 1% for PID. Using sev¬ 
eral large inclusive kaon and 73-decay samples, we find 
a systematic uncertainty for A c h of 1.1%, due mainly to 
the dependence of reconstruction efficiency on the charge, 
for the high momentum pion from 73 + —> t]tt + , pK + and 
rf 7r + . The corresponding number for the softer charged 
pion from the p + in B + —> r)p + is 2%. 

In this Letter, we have presented improved mea¬ 
surements of branching fractions for six charmless 73- 
meson decays. All branching fractions are in agreement 
with theoretical predictions. The previously unobserved 
B+ —» r]p + and 73 + —> ?/7r + decay modes are seen with 
significance 4.7er and 5.4cr, respectively. For the charged 
modes, we also determine the charge asymmetries. These 
are found to be consistent with zero within their uncer¬ 
tainties. 
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